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AMEND3MENTS TO THE CLAIMS 
Please amend the claims as indicated. 



1 1 . (withdrawn) A method of modeling faulting and fracturing in a subsurface 

2 volume of the earth comprising: 

3 (a) selecting a mode of definition of a subsurface model, said mode of 

4 definition selected from (i) an aerial mode wherein the model comprises a 

5 plurality of nodes in a horizontal plane interconnected to each other and to 

6 a substrate, (ii) a cross sectional mode wherein the model comprises a 

7 plurality of nodes in a vertical plane interconnected to each other and a 

8 substrate defining the edges of the model, and, (iii) a 3-D mode wherein 

9 the model comprises a plurality of nodes interconnected to each other and 

10 to a substrate defining the edges of the model; 

11 (b) defining said subsurface model including specifying material rock 

1 2 properties within the subsurface volume; 

] 3 (c) specifying an initial deformation pattern; and 

1 4 (d) using a dynamic range relaxation algorithm to find a force equilibrium 

15 solution for said subsurface model and said initial deformation pattern 

16 giving a resulting deformed model including fracturing. 
17 

1 2. (withdrawn) The method of claim 1 > wherein selecting said mode of definition, 
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2 defining a subsurface model, and specifying said initial deformation pattern 

3 further comprises using a graphical user interface. 
4 

1 3. (withdrawn) The method of claim 1 , wherein said nodes are arranged in a grid 

2 that is one of (i) a triangular grid, and, (ii) a random grid. 

3 

1 4. (withdrawn) The method of claim 1 , wherein said nodes are interconnected by 

2 springs, and defining said subsurface model further comprises defining a normal 

3 force associated with each spring. 
4 

1 5. (withdrawn) The method of claim 4, wherein defining said subsurface model 

2 further comprises a substrate attachment force associated with each node that is 

3 attached to said substrate. 
4 

16. (withdrawn) The method of claim 1 , wherein said nodes are interconnected by 

2 beams, and defining said subsurface model further comprises defining a normal 

3 force and a shear force associated with each beam . 
4 

1 7. (withdrawn) The method of claim 1 , wherein specifying said initial deformation 

2 pattern further comprises performing a reconstruction based at least in part upon 

3 an observed large-scale deformation corresponding to said subsurface volume. 
4 
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1 8. (withdrawn) The method of claim 7, wherein said reconstruction is a palinspastic 

2 reconstruction. 
3 

1 9, (withdrawn) The method of claim 7, wherein obtaining said initial deformation 

2 pattern further comprises: 

3 (i) obtaining a trial deformation pattern from said observed large scale 

4 deformations, 

5 (ii) applying an anticipate method to said model using said trial deformation, 

6 giving an approximate deformation result wherein said approximate 

7 deformation result is exclusive of fractures or faults; and 

8 (iii) updating said trial deformation based on a comparison of said approximate 

9 deformation result and said observed large scale deformation thereby 
10 giving said initial deformation pattern, 

11 

1 1 0. (withdrawn) The method of claim 7 further comprising conditioning said 

2 subsurface model thereby increasing the likelihood of said resulting deformed 

3 model including said observed large scale deformations, said conditioning 

4 including a weakening of bonds between adjacent ones of saia plurality of nodes 

5 over at least a portion of the subsurface model. 
6 

1 11. (withdrawn) The method of claim 1 , wherein using the dynamic range relaxation 

2 algorithm further comprises applying said initial deformation model to said 
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3 substrates in a plurality of steps, each step comprising a applying specified 

4 fraction of the initial deformation to said substrates and determining if any bonds 

5 between the nodes have been deformed beyond a breaking point and identifying a 

6 subset of the bonds that have been so deformed. 
7 

1 12. (withdrawn) The method of claim 1 1 , wherein applying the dynamic range 

2 relaxation algorithm further comprises iterativeiy breaking the one bond of the 

3 subset of bonds that has been deformed the most and applying a relaxation 

4 algorithm to the remaining unbroken bonds. 
5 

1 1 3- (withdrawn) A method of simulating faulting and fracturing in a subsurface 

2 volume modeled by a plurality of interconnected nodes due to an initial 

3 deformation pattern applied to boundaries of said subsurface volume, the method 

4 comprising the following steps: 

5 (a) applying a fraction of said initial deformation to said boundaries; 

6 (b) using a dynamic range relaxation algorithm (DRRA) to find a force 

7 equilibrium solution for said applied fractional initial deformation and 
g identifying bonds between said plurality of interconnected nodes 

9 susceptible to breakage, wherein identifying bonds susceptible to breakage 

1 0 further comprises comparing a deformation of each bond in said force 

I ! equilibrium solution to at least one predetermined breakage threshold 

1 2 associated with each of said bonds; 
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1 3 (c) if in step (b) no bonds susceptible to breakage are identified, increasing 

14 said initial fractional deformation of said boundaries and increasing said 

1 5 fraction of said initial deformation and iteratively repeating steps (a) - (c ) 

16 until said fraction equals one; and 

17 (d) if in step (b), at least one bond is susceptible to breakage is identified, then 

18 breaking the one of the identified bonds susceptible to breakage whose 

19 deformation exceeds its at least one associated breakage threshold the 

20 most, and repeating step (b) 
21 

1 14. (withdrawn) The method of claim 1 3, wherein using a DRRA further comprises: 

2 (i) relaxing the plurality of interconnected nodes according to a single over- 

3 relaxation step to give a relaxed position of said plurality of nodes; 

4 (ii) identifying a first subset of the plurality of nodes that move further than a 

5 relaxation threshold; 

6 (iii) if said first subset of the plurality of nodes is empty, using said relaxed 

7 positions as said force equilibrium solution. 
8 

1 15. (withdrawn) The method of claim 1 4, wherein said first subset of nodes at step 

2 (iii) is not empty, and using the DRRA further comprises: 

3 I. relaxing sequentially each of the nodes in said first subset of nodes and 

4 identifying a second subset of the plurality of nodes comprising those of 

5 the first subset of nodes, and each node connected thereto, that move more 
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6 than a relaxation threshold ; 

7 H. after step II, interchanging the nodes in the first and second subset of 

8 nodes; and 

9 III. interatively repeating steps 1 and II until the first subset of nodes is empty. 
10 

1 1 6. (withdrawn) The method of claim 13, wherein the connections between the 

2 plurality of interconnected nodes comprise a plurality of springs and at least one 

3 associated breaking threshold is an extensional breaking threshold. 
4 

1 1 7. (withdrawn) The method of claim 1 6, wherein the plurality of extensional 

2 breaking thresholds comprise one of (i) a Gaussian distribution characterized by a 

3 mean value and a standard deviation, and, (ii) a Weybull distribution. 
4 

1 18. (withdrawn) The method of claim 1 3, wherein the connections between the 

2 plurality of interconnected nodes comprise a plurality of beams and at least one 

3 associated breaking threshold comprises an extensional breaking threshold and a 

4 shear breaking threshold 
5 

1 19. (withdrawn) The method of claim 1 8, wherein the plurality of extensional 

2 breaking thresholds comprise one of (i) a Gaussian distribution characterized by a 

3 mean value and a standard deviation, and (ii) a Weybull distribution, and the 

4 plurality of shear breaking thresholds comprise one of a Gaussian distribution and 
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5 a Weybull distribution. 
6 

1 20. (withdrawn) The method of claim 13, wherein the subsurface volume further 

2 comprises a plurality of regions, each of said plurality of regions characterized by 

3 an associated materials having material properties. 
4 

1 21 . (withdrawn) The method of claim 20, wherein said associated materials are 

2 selected from the group consisting of (i) salt, and, (ii) a rock. 
3 

1 22. (withdrawn) The method of claim 1 3, wherein the plurality of interconnected 

2 nodes constitute an aerial network and said connection between said 



3 interconnected nodes is selected from the group consisting of (i) springs, and, (ii) 

4 beams, and wherein the boundaries further comprise a plurality of substrate 

5 nodes, said plurality of substrate nodes attached to proximate nodes of the aerial 

6 network by the same type of connection as the connection between the 

7 interconnected nodes. 
8 

(withdrawn) The method of claim 13, wherein the plurality of interconnected 
nodes constitutes a 2-D cross section and wherein the boundaries comprise a 
plurality of discs, each of said plurality of discs experiencing at least one of (i) an 
attractive force towards, and, (ii) a repulsive force away from at least one of the 
plurality of interconnected nodes. 
6 
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1 24 , (withdrawn) The method of claim 1 3 , wherein the plurality of interconnected 

2 nodes constitutes a 3-D network and wherein the boundaries comprise a plurality 

3 of spheres, each of said plurality of spheres experiencing at least one of (i) an 

4 attractive force towards, and, (ii) a repulsive force away from at least one of the 

5 plurality of interconnected nodes. 
6 

1 25. (withdrawn) The method of claim 23 further comprising checking a distance 

2 between pairs of said plurality of discs to a predetermined threshold after step (a) 

3 of claim 1 3 and adding additional discs to said boundaries if said distance exceeds 

4 said predetermined threshold. 

5 

1 26. (withdrawn) The method of claim 24 further comprising checking a distance 

2 between pairs of said plurality of spheres to a predetermined threshold after step 

3 (a) of claim 13 and adding additional spheres to said boundaries if said distance 

4 exceeds said predetermined threshold. 
5 

1 27. (withdrawn) The method of claim 1 8 further comprising applying a von Mises 

2 failure criterion. 
3 

1 28. (withdrawn) The method of claim 1 3 further comprising preconditioning the 

2 model to increase the likelihood of fracturing of the bonds between the plurality 

3 of interconnected nodes in the vicinity of specified locations. 



4 
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1 29. (withdrawn) The method of claim 28, wherein said preconditioning further 

2 comprises reducing the predetermined breakage threshold of those of said bonds 

3 in the vicinity of said specified locations. 
4 

1 30, (withdrawn) The method of claim 29 3 wherein said specified locations are 

2 obtained by a geologic reconstruction based on observed large scale deformations. 
3 

1 3 L (withdrawn) The method of claim 29, wherein said specified locations further 

2 comprise one of (i) piecewise linear curves, and (ii) piecewise linear surfaces. 
3 

1 32. (withdrawn) The method of claim 29, wherein reducing the strength of said bonds 

2 in the vicinity of the specified locations is based at least in part on a linear 

3 function of distance from said specified locations. 
4 

1 33 . (withdrawn) A method of simulating deformation without faulting and fracturing 

2 due to an initial deformation pattern applied to boundaries of a subsurface volume 

3 modeled by a plurality of interconnected nodes, the method comprising: 

4 (a) defining a plurality of boundary nodes on a boundary of said subsurface 

5 volume wherein said initial deformation pattern is applied; 

6 (b) defining an initial and a final position for each of said plurality of 

7 boundary nodes and a displacement there between; 

8 (c) determining a distance from each of the plurality of interconnected nodes 

9 to the final positions of the plurality of boundary nodes; 
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10 (d) determining a displacement for each of the plurality of interconnected 

1 1 nodes as a combination of said displacement of said boundary nodes 

12 weighted by a weighting function related to said distance from each of the 

13 plurality of interconnected nodes and the final positions of the boundary 

14 nodes. 
15 

1 34, (withdrawn) The method of claim 33, wherein said weighting function includes an 

2 exponential factor related to said distance fiom each of the plurality of 

3 interconnected nodes and the final positions of the boundary nodes. 
4 

1 35. (withdrawn) The method of claim 34, wherein said weighting function further 

2 includes a matrix whose coefficients are obtained by solving a plurality of 

3 equations including said displacement of the boundary nodes wherein the plurality 

4 of equations is three times the plurality of said boundary nodes. 
5 

1 36, (currently amended) A graphical user interface (GUI) for displaying and 

2 manipulating a model of interconnected nodes for simulating fracturing and 

3 faulting in a subsurface volume of the earth, comprising: 

4 (a) a first module for pr e s e nting which presents in a portion of a computer 

5 screen, a first graphical image representative of a plurality of 

6 interconnected nodes of the model; 

7 (b) a second module for d e fining which defines m aterial properties of the 
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8 model defined by the plurality of interconnected nodes; 

9 (c) a third module for defining which defines a n initial deformation pattern 

10 applied to boundaries or the substrate of said plurality of interconnected 

11 nodes; 

12 (d) a fourth module for d e fining which defines p arameters of a simulation 

1 3 process including a dynamic range relaxation algorithm for simulating a 

14 response of said model to said initial deformation pattern. 
15 

1 37. (currently amended) The GUI of claim 36, wherein said first module is capable 

2 of presenting displays said first graphical image at least one of (i) a planar view, 

3 (ii) a cross-sectional view, (iii) a 2D view, and (iv) a 3D view. 
4 

1 38. (original) The GUI of claim 36 further comprising the input of a random number 

2 seed for the random number generator used for setting up at least one of (i) a 

3 geometry of said interconnected nodes in said model, and (ii) breaking thresholds 

4 associated with links between pairs of interconnected nodes in said model. 
5 

1 39. (currently amended) The GUI of claim 36, wherein said fourth module further 

2 comprises an editor for setting which defines at least one of (i) a relaxation 

3 threshold for said simulation, (ii) an over-relaxation factor for said simulation, 

4 (iii) a maximum movement during said simulation, (iv) a time step for said 

5 simulation, (v) an angular relaxation factor for said simulation, and (vi) an angular 

6 over-relaxation factor for said simulation. 
09/923>048 14 
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7 

1 40. (currently amended) The GUI of claim 36, wherein said third module further 

2 comprises an editor for defining which defines a deformation that is at least one of 

3 (i) a localized extension, (ii) uniform extension, (iii) uniform compression, (iv) a 

4 uniform right lateral shear, (v) a uniform left lateral shear, (vi) rotation, (vii) a 

5 deformation region in areal simulation mode or on the lowest plane in 3-D 

6 simulation, (viii) translation to a deformation region, and (ix) a rotation to a 

7 deformation region. 
8 

1 4 L (currently amended) The GUI of claim 36, wherein said second module further 

2 comprises a material editor for dofining which defines a t least one region selected 

3 from the group consisting of (i) a rock region in said model, and (ii) a salt region 

4 in said model 
5 

1 42. (currently amended) The GUI of claim 4 1 , wherein the material editor further 

2 eomprisofl defining defince, for each link associated with each pair of the plurality 

3 of interconnected nodes for the at least one region, properties a property selected 

4 from (A) an extensional breaking threshold for said link, (B) a shear breaking 

5 threshold for said link, (C) a linear force constant for said link, and (D) a shear 

6 force constant for said link. 
7 

1 43. (currently amended) The GUI of claim 36, wherein said first module further 

2 comprises an editor for controlling which controls t he display of at least one of (i) 
09/923,04$ 15 



PAGE 16(25 * RCVD AT 7/25/2005 12:28:04 PM [Eastern Daylight Time] 1 SVfcUSPTO-EFXRF-6/29 1 DNS:27383QQ 1 CSID:7132668510* DURATION (mm-ss):05-38 



MfiDftNMOSSMflNSRIRfiM Fax : 71 32668510 



Jul 25 '05 11:34 



P. 1? 



3 faulting resulting from said simulation process, and (ii) a stress s fresses resulting 

4 from said simulation process. 
5 

1 44. canceled 

2 

1 45 . (currently amended) The GUI of olaim - H claim 36. further comprising a module 

2 to display which displays a g raphical images image o f the model output after 

3 simulating deformation. 



4 

1 46. canceled 

2 

1 47, canceled 

2 

1 48. canceled 

2 

3 49. canceled 
4 

5 50. canceled 

6 

1 51. (currently amended) The GUI of okumSO claim 40 further comprising an editor 

2 which, for defining deformation, restricts the motion of the cursor to the 

3 coordinate system of the model. 

4 
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1 52. canceled 

2 

1 S3. canceled 

2 

1 54. canceled 
2 



1 55. (currently amended)The GUI of claim 44 claim 36. wherein a module is-capable 

2 of controlling for controls d isplay said ofa^raphical images image o f a model 

3 that is gridded. 
4 

1 56, (previously presented) The GUI of claim 55, wherein a modulo io capable of 

2 controlling for display said graphical imag e s of a - die model &at is gridded as on 

3 at least one of (i) a triangular grid, (ii) a tetrahedral grid, (iii) a rectangular grid, 

4 and (iv) a random grid. 
5 

1 57. (currently amended) Hie GUI of claim 11 clam36> wherein a modulo io capabl e 

2 of controlling for display said graphical imagoo of a model with thejiodes are 

3 interconnected by at least one of (i) springs, and, (ii) beams. 

4 

1 58. (currently amended) The GUI of efekjHM- claim 36, wherein a modulo is capable 

2 of controlling for display said graphical images of a mod e l with interconnected 

3 the nodes having frave associated forces that are at least one of (i) normal forces, 

4 (ii) shear forces, (iii) attractive forces, (iv) repulsive forces, and, (v) substrate 
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5 attachment forces. 
6 

1 59, (currently amended) The GUI ofe k*im 44 claim 36, wherein a module is c ap ab le 

2 of controlling for display an displays the initial deformation pattern. 
3 

1 60. (currently amended) The GUI of claim 44 claim 36. furth e r comprising a module 

2 wherein extensional breaking thresholds and shear breaking thresholds between 

3 the i nterconnected nodes may be defined as statistical distributions. 
4 

1 61 . (currently amended) The GUI of oloim H claim 36, furth e r comprising a modul e 

2 wherein breaking thresholds between interconnected the nodes mayb e are defined 

3 as bygo ne of (i) a Gaussian distribution characterized by a mean value and a 

4 standard deviation, and (ii) a Weybull distributio n, and wher e in t he plurality of 

5 shear breaking thresholds compris e one of a Gaussian distribution and a Woybull 

6 distribution . 
7 

1 62. (currently amended) The GUI of olaim - 4 4 claim 36 further comprising a module 

2 for cheoking which compares a distance between pairs a pair of said 

3 interconnected nodes to a predetermined threshold and adding adds an additional 

4 a odcs nodej o amodel boundari e s boundary if said distance e xceeds a exceeds the 

5 predetermined threshold. 
6 

1 63. (currently amended) The GUI of claim 4 4 claim 36 w herein a module is capabl e 
09/923,048 18 
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2 of d e fining modol preconditioning preconditions the model t o increase &e a 

3 likelihood of fracturing of fte bonds between a plurality of interconnected nodes 

4 in a vicinity of specified locations. 
5 

1 64, (currently amended) The GUI of claim 63 wherein said preconditioning further 

2 comprises reducing fee predetermined breakage threshold of those of said bonds 

3 in the vicinity of said specified locations. 
4 

1 65. (previously presented) The GUI of claim 63 wherein said specified locations 

2 further comprise one of (i) piecewise linear curves, and (ii) piecewise linear 

3 surfaces. 
4 

1 66. canceled 

2 

1 67. (currently amended) The GUI of claim 44 claim 36, further comprising a module 

2 to di s pl ay which displays f ault surfaces projected on to at least one of (i) a 

3 substantially planar view, (ii) a cross-sectional view, (iii) a 2D view^ and (iv) a 3- 

4 D view* 
5 

1 68. (currently amended) The GUI of claim 41 claim 36. further comprising a module 

2 to disp l a y which displays g raphical images wherein the interconnected nodes are 

3 displayed using a rendering quality dependent upon a predetermined deformation 

4 size and degree. 
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5 

1 69. (currently amended) The GUI of claim 44 claim 36. further comprising a module 

2 to display which displays g raphical images wherein the interconnected nodes are 

3 displayed dependent upon whether a predetermined deformation threshold has 

4 been reached. 
5 

1 70. (previously presented) The GUI of claim 69 wherein said predetermined 

2 deformation threshold is dependent upon at least one of (i) a specified number of 

3 broken connections between nodes and (ii) a specified change in length between 

4 nodes. 
5 

1 71. (currently amended) The GUI of claim 11 claim 36, further comprising a module 

2 to display which displays g raphical images of tike astress distribution within a 

3 deformed network of said interconnected nodes. 
4 

1 72, (previously presented) The GUI of claim 7 1 wherein said stress distribution is 

2 displayed using a color bar and scale, the color dependent on at least one of (i) a 

3 scalar stress quantity for each node, (ii) a scalar stress quantity for each node 

4 interconnection, (iii) a scalar stress quantity for a grid^ and (iv) a scalar stress 

5 quantity averaged over at least two grids. 
6 

1 73 . (currently amended) The GUI of claim 9 claim 36, further comprising a module te 

2 display which displays g raphical images of deformation as an animation series. 
09/923,048 20 
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3 

1 74. (currently amended) The GUI of claim 44 claim 36. further comprising a module 

2 to display which displays imported graphical images superimposed with said 

3 model of a plurality of interconnected nodes for simulating to illustrate f racturing 

4 and faulting in a subsurface volume of the earth. 
5 

1 75. (previously presented) The GUI of claim 74 wherein said imported images aaay 

2 fee^ne are selected from the group consisting o f (i) a seismic section, (ii) a 

3 geologic cross section and (iii) an arbitrary earth model. 
4 
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